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Thanks to dimensional reduction, the contributions to the hot QCD pressure coming from so- 
called soft modes can be studied via an effective three-dimensional theory named Electrostatic 
QCD (spatial Yang-Mills fields plus an adjoint Higgs scalar). The poor convergence of the per- 
turbative series within EQCD suggests to perform lattice measurements of some of the associated 
gluon condensates. These turn out, however, to be plagued by large discretization artifacts. We 
discuss how Numerical Stochastic Perturbation Theory can be exploited to determine the full lat- 
tice spacing dependence of one of these condensates up to 4-loop order, and sharpen our tools on 
a concrete 2-loop example. 



The XXV International Symposium on Lattice Field Theory 
July 30-4 August 2007 
Regensburg, Germany 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Towards 4-loop NSPT result for a 3d condensate-contribution to hot QCD pressure 



C. Torrero 



1. Motivation 

The pressure p(T) of hot QCD (with a temperature T S> 150 MeV) is of crucial importance in 
different contexts. For instance, in cosmology, the cooling rate of the Universe is given by 



1 dr V247T y/e{f)s(T) 



(1.1) 



T dt m P i c(T) 

where T is the temperature, t the age of the Universe, m P1 the Planck mass, and 

*(T)=y(r), e{T) = Ts(T)-p{T) , c = e'(T). (1.2) 

On the other hand, in /jeavj ion collisions, the system evolves according to 

T^ v ^[p(T)+e(T)]u^u v -p(T)rj^ v , 9^ = 0, (1.3) 

where 7^ v is the energy-momentum tensor, r]^ v the Minkowski metric, and the flow velocity. 
Finally, from a theoretical point of view, the pressure is proportional to the number of effective 
degrees of freedom, and is therefore a good observable to characterize hot QCD matter. 



2. Theoretical setup — part I (continuum) 

In general, the determination of the QCD pressure is a non-trivial task. In spite of our re- 
striction to the deconfined phase, T 3> 150 MeV, where perturbation theory should in principle be 
applicable, it is of limited use in practice, because of the very slow convergence of the perturba- 
tive series [Q]. At the same time, first principles lattice simulations are also difficult at temperatures 
above about T ~ 1 GeV because the system then develops a scale hierarchy, g 2 T jn <gT < nT , 
where g is the QCD gauge coupling (for an idea to possibly overcome this limitation, see ref. 

A suitable strategy to tackle the computation of p(T) (and of many other observables) in this 
situation is given by Dimensional Reduction [|], ||, gj: it consists of integrating out "hard modes", 
with momenta k ~ nT, from the four-dimensional (4d) QCD, to arrive at an effective description in 
terms of so-called Electrostatic QCD (EQCD) i.e. a three-dimensional (3d) Yang-Mills theory 
plus an adjoint Higgs field Aq(jc). The action of EQCD is given by 

S EQCD = / d 3 * j^Tr[^(x)] +Tr[D ; ,A (x)] 2 + m 2 E Tr[A 2 (x)} + ^(Tr[A 2 (x)} ) 2 | , (2.1) 

where is the 3d field strength tensor; D,- the covariant derivative; Ao = L| =1 Aq T b with 
Tr [T A T B ] = \8 AB ; and we implicitly assume the use of dimensional regularization. The tempera- 
ture T now enters only via the parameters m E , A E and g E , where g E is the EQCD gauge coupling. 

It turns out that perturbation theory within EQCD might converge very slowly [Q] (see, how- 
ever, ref. [§]). One can then switch to a numerical measurement of the partial derivatives of the 
pressure with respect to the parameters of the EQCD action, the so-called condensates: after sub- 
tracting the proper counterterms, the results are extrapolated to the continuum and then numerically 
integrated to finally get p(T) It turns out that carrying out the continuum extrapolation is dif- 
ficult because of large 0(a) discretization effects: thanks to super-renormalizability, these effects 
are however purely perturbative in nature, and our aim is to determine them up to 4-loop level. 
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3. Theoretical setup — part II (lattice) 

The EQCD action on the lattice may be written as 

S m = P £ (l - 1 Re Tr [/> 7 (*)])- 

X, i<j 

- 2 £ Tr [0(x) Ui(x) 4>(x+i) U?(x)] + 

X, ( < j 

+ £|a(j3,A,am E )Tr^ 2 (x)] + A(Tr [0 2 (jc)]) 2 | , (3.1) 

where j8 = 6/ag 2 , Pjj is the plaquette, t/,- is the link variable, = ^/oAq, X = aX E , and [ |To| ] 

, f (am E ) 2 / 5, \ 3.175911525625 
a(P,X,am E ) = 6|l + ^-- (6 + -A/3J ^ 



87T 2 j3 : 



WXP-^X 2 pA f In j8 +0.08849^) + + 36.130 | 



(3.2) 



The quantity under inspection is the derivative of p(T) with respect to [am E ) 2 \ this yields (Tr [0 2 ] ), 
which can be expanded as 



1 1 A 

( Tr [ (j) 2 } ) = doo + dm ^ + d n X + d 2 o + rf 2i j + d%%X 2 + 

I X X 2 ( X" 

+ d 30 p+d 3i -^ + d 32 j+d 33 X 3 + O^j^ 



(3.3) 



The coefficients doo, d\o, d\\, d 2 \ and d 22 are known analytically, for instance (N = lattice extent) 



d u = 40 
1 

X Jp 



3.175911525625 1 



An 

N-l N-l N-l 

III 

k l= k 2 =0 k 3 =0 



N-l N-l N-l 
III 

Jl=0 H2- 

1 



1 



A^ 3 n % n % n % 4L 3 =1 sin 2 (^) + (am E ) 2 

~i 2 



4L 3 =1 sin 2 (f ) + Kf 



(3.4) 



but the others have to be determined: those with the biggest impact are expected to be the 3-loop 
and 4-loop coefficients independent of X, i.e. d 2 o and d 3 o, respectively. 



4. Numerical setup 

The perturbative study is concretely carried out by means of Numerical Stochastic Perturba- 



tion Theory (NSPT) [11 ]. (Incidentally, it would also be interesting to pursue the same computation 



with standard techniques [12].) Its origins lie in Stochastic Quantization [13], based on introducing 
an extra coordinate t and an evolution equation of the Langevin type, namely 



d t (j)(x,t) = -<fyS[0] + rj(x,t) , 



(4-1) 
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Figure 1: Analytical and numerical results for the coefficient d\ \ vs. N at am E = 0.4. The results agree with 
each other within statistical errors. 



where f]{x,t) is a Gaussian noise. The usual Feynman-Gibbs path integral is recovered by averag- 
ing over the stochastic time, 

Z- l J [D<l>}0[Hx)]e- Smx)] =1^7 jfV<0fo,(x,O]> T , • (4-2) 
NSPT can now be introduced by expanding the variables as 

<Hx,o-^£^ w Cm), (4-3) 

k 

where go is some small coupling. This results in a hierarchical system of differential equations that 
can be numerically integrated by discretizing the stochastic time, as t = nx, where X is a time step. 

A similar construction holds also for the gauge degrees of freedom Uj(x), for which the 
Langevin equation reads 

d t U v =-i(vS[U ri ] + T 1 ^Ur 1 , (4.4) 

in order to assure a correct evolution within the group. The perturbative expansion is then a double 
expansion in j3 and A, to obtain the previously-written series of (Tr [(f) 2 ] }. 

In practice, every variable evolves according to its Langevin dynamics for different values of 
t; a measurement of Tr [<p 2 ] is performed at every time step (once thermalization has been reached); 
and finally one extrapolates to T = (this last step is necessary since the correct probabilistic weight 
at equilibrium is recovered only in the limit z — > 0). This procedure is then repeated after changing 
the parameters of the action. 



5. Preliminary results 

Our approach involves three different extrapolations / interpolations in total: first, the above- 
mentioned extrapolation to x — >• 0; second, an extrapolation to infinite volume (N — > °°); third, an 



4 



Towards 4-loop NSPT result for a 3d condensate-contribution to hot QCD pressure 



C. Torrero 



= 0.3 



1.0 



0.20 



0.18 



0.16 



0.14 



0.00 



■ 


exact 


o 


NSPT 



0.05 0.10 
IfN 



0.15 



0.129 



0.128 



0.127 



0.126 



■ 


exact 


o 


NSPT 



0.00 



0.05 0.10 
1/N 



0.15 



Figure 2: Results for the coefficient d\\ vs. l/N at am E = 0.3 and 1 .0. Note the very different resolutions of 
the vertical axes. 

interpolation between the different am E simulated, in order to obtain the desired 4-loop predictions 
at any finite am E . 

As a check of the first of these extrapolations, we compare the numerical outputs for d\\ with 
its known values (eq. (^4j)) at fixed am E and lattice extent N: this is done in Fig. 1. 

The next task is to extrapolate to infinite volume. As Fig. 1 shows, finite-volume effects 
become exponentially small at large volumes. However, the volume required for this grows as the 
mass am E decreases. This is illustrated in Fig. 2 for two masses, one larger than in Fig. 1, and 
one smaller. For the large mass, a plateau can be reached allowing for a reliable infinite-volume 
extrapolation, while for the smaller mass, the largest volumes we can afford are not yet large enough 
to reach a plateau. 

In order to deal with this situation, we adopt the following approach. Let us consider a mass 
like am E = 0.8: numerical evidence shows that this one has a reasonable plateau at affordable 
N for all the coefficients djj, but still the behaviour of the data is not too flat (i.e., some volume 
dependence is detectable). One can then extract an infinite- volume value dij(°°) by fitting a constant 
to data in the range of the plateau, and subtract it from the data in order to obtain the quantities 
(L = Na) 



8ij{m E L) 



: dij{m E L) 



■dijH 



(5.1) 



Subsequently, one can try to obtain a reasonable interpolating fit fij(m E L) for g, ; -(m E L), allowing 
to go also to other values of m E L than those simulated at am E = 0.8. After this, one can go back to 
the other masses m' E , and take a finite-size scaling ansatz of the form 



dij(m' E L) 



: dij(oo) +Aij{am' E ) x fij{m' E L) , 



(5.2) 



where dij(m' E L) are the direct measurements at the mass m' E , and dij(oo) and Ajj(am' E ) are volume- 
independent fit coefficients. Test results for d\\ (°°) obtained this way are shown in Fig. 3. 
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Figure 3: Comparison between exact values (dashed curve) and numerical infinite-volume extrapolations 
(open symbols) as a function ofam E . The continuum value at am E — has been extracted from ref. The 
dotted line shows a fit through the open symbols, constrained to go through the continuum point ( which is 
known for all dij). The agreement between the exact and fitted curves is satisfactory. 

6. Conclusions 



Preliminary studies of the 2-loop coefficient d\\, where analytical values are also available, 
show that our general approach works. The next goal is to finalise the analysis for the most impor- 



tant 3-loop coefficient J20 and 4-loop coefficient J30 Q15[]. The quality of our data is good, so we 
expect to be able to extract these with small errors in the same mass range as in Fig. 3. This should 
allow to re-analyse the Monte Carlo data of ref. ^ with significantly reduced systematic errors 
from the continuum extrapolation. Combining with refs. fli6[ [l7[ , |18| ], all "soft" contributions to the 
hot QCD pressure would then be under reasonable control. At the same time, the determination of 
the 4-loop "hard" contributions remains an open challenge; toy model computations in scalar field 



theory have suggested, however, that it can be tackled with some effort [19]. 
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